On the Raman shift in nanosized crystals 
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Abstract. An analytical form of the Raman shift dependence on size of nanocrys- 
tals is presented. Based on the hard confinement model, it works in those cases 
where the average phonon curve shows a quadratic dependence on the phonon 
quasi- momentum in the range of interest. 
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Raman spectroscopy is an important tool in gathering information about molec- 
ular and crystal vibrational properties. Raman spectra can show significant 
changes when sizes of investigated crystals go down to the nanometer scale. 
Strain, stress and non-stoichiometry can be invoked to explain these changes 
[1,2]. However, many published works have shown that in most cases the main 
cause of spectra changes is to be ascribed to the phonon confinement (PC) ef- 
fect [3 — 9] . Actually, the confinement of the go = phonon wave function in 
nanometric-sized crystals makes accessible to Raman investigation a significant 
portion of the Brillouin Zone (BZ), whose extension increases as the crystal size 
decreases [3 — 5]. Besides some specific cases [5, 10, 11], a direct relation, even 
when approximated, which explicitly shows the connection between Raman shift 
and crystal size, depending on suitable model parameters, is still not available. 
In this letter we present a detailed numerical analysis of the hard confinement 
model (HC) which yields an unexpectedly simple analytical form capable of 
closely reproducing the model predictions. It uses the quadratic approximation 
of the crystal phonon curve ui{q) which, in several practical cases, holds down 
to a few nanometers of crystal sizes. 

Formally, the PC effect can be explained by considering the modulation of 
the phonon wave function in the infinite crystal with a suitable "weighting" 
function [4]. On expansion by Fourier integrals, we can calculate the first order 
Raman spectrum of a nano-sized crystal as 

J(w) oc [n(w) + 1] / dq (1) 
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where C(q) stand for Fourier coefficients of the go = phonon wave function, 
uj(q) for the phonon dispersion curve, Tq for the intrinsic Raman band linewidth 
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and [n(oj) + 1] for the Bose-Einstein factor [4]; the latter being usually disre- 
garded for lineshape characterization or nanoparticlc sizing. Some simplifica- 
tions are required in order to handle eq (P). A drastic but commonly accepted 
assumption considers an isotropic dispersion in a spherical BZ [4]. In this case, 
the function uj(q) represents an average dispersion curve [6]. As concerns phonon 
confinement , the HC model assumes a gaussian weighting function for spher- 
ical nanocrystals. Accordingly, by disregarding unessential factors, eq (fTJ) can 
be rewritten as [6, 8, 15] 

/(a,) = T^i-Wfi«« a )<? dq (2) 

J [c-^)] 2 + (r /2)2 

where L stands for the particle size and a for the lattice parameter. The latter 
equation marks the start line of the present work. 

If the quadratic approximation holds the phonon curve u>(q) can be repre- 
sented as 

uj(q) =uj a + Aq 2 a 2 /2 (3) 

where A stands for a suitable curve parameter. Equation |3]) provides an univo- 
cal correspondence between peak frequency 2 and an effective quasi-momentum 
q which allows for u = u>(q). Thus, we should search for an expression of w by 
investigating the proper definition of q. 

For convenience , let us define the function 

G(q, L) = q 2 exp(-q 2 L 2 /16 vr 2 ) (4) 

7T / a 

so that, I(uj) = J G(q,L)/{[^~u;(q)} 2 + (T Q /2) 2 }dq . By taking into ac- 
o 

count that at the maximum of I(u>) the condition [d I{w)/dw\ ~ = holds, 
we obtain the exact equation 

G^ L )p- u m 2dq=0 (5) 



{[5-c(g)] 2 + (r /2)2}^ 



Of course, in the framework of the HC model the latter equation is general, that 
is, not related to a special law for to(q). Thus, for a given Tq and w(q), we can 
define the distribution-like function 

G(q,L) . 7° G(q,L) 



MQ,L) = - f-t -2/ /- f-^ -sdq (6) 

{p-c(g)] 2 + (r /2)2} { {p-c( g )] 2 + (r /2)2} 



so that 

7r/a 



Zj= \ Fz(q,L)u(q)dq (7) 
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Now, by inserting form ([3]) into eq. |(7J) we can make explicit the meaning of the 
effective quasi-momentum q , that is, 



LO = OJQ + 



■k I a 

A f A 

—a 2 / F%(q, L)q 2 dq = lu q + ^q 2 



(8) 



The analytical resolution of cq. ([5]) is quite concealed and its reduction 
to a differential equation, even if simplified, is not helpful for the purpose. 
However, for the sake of discussion, it is worth while to handle eq. ((5} in a more 
explicit form. By taking into account that G(n/a,L) k ( L > a) and after 
integrations by parts, derivations with respect to L and by disregarding small 
terms we obtain 



em 



A(i) 



9 [1 - &{L)Y + A(L) 



d@(L) 



dL 



where 



and 



O(L) = 1 



7T J qa 

A(L)= J F s {t,L) 
o 



24tt 2 



(l-C 2 )/e 2 



{(l-e 2 ) 2 + (Wax)4} 



0) 
(10) 
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with x = {\A\ /ro) 1 ^ 2 . Due to the gaussian factor present in F%(£, L) , if L is 
sufficiently large with respect to a , the upper integration limit can be considered 
as infinite so that A could be conveniently presented as a function on 77 = L jxa 
rather that on L. Thus, to avoid confusion, further on we use special symbols 
©s( 7 ?) = an d As(?7) = A(L) even when small size are considered. 

It is evident from eq. that as L increases As (77) decreases because q —> 0. 
Thus, for large L the effective quasi-momentum decreases as q 2 w 24ir 2 /L 2 and 
As(rj), as well as ®s(il)> decreases as ~ l/r/ 4 . In the range of small sizes, 
eq. @ cannot be used. However, from eq. ([5]) we expect that for vanishing 
crystal sizes Qsiv) ~ * 1 because of q < n/a. In the limiting case of q — > n/a 
we expect 65(77) ~ exp(— x 2 r/ 2 /2A). A more detailed description of 85(77) is 
obtained from the numerical calculation of q. Results are shown in the inset of 
Fig.l as obtained for different x values. For comparison, segments of parabola 
approaching the curves —In [65(77)] are shown as well. Precisely, parabola have 
equations f(rj) = mx 2 rj 2 /24 where m = 0.65, m = 0.8, m = 1 and m = 1 for 
x = 0.5, x = 1, x = 2 and x — 3, respectively. 

Usually A > To so that our interest will be focused on the x > 1 curves. For 
77 > 7 the latter curves are independent of x and quickly approach the straight 
line g(j]) = 77/ 12, maintaining this trend up to 77 w 60 where 65(77) << 1 
(Fig.l). On account of eq. (flT)|) , we should not be concerned about the deviation 
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of — In [0,5(77)] from the linear law since for 77 > 60 we have q 2 24ir 2 /L 2 . Thus 
we can be confident of the approximation 

f ap = [1 - exp(-L/12az)] (12) 

and, consequently, 

Aa 2 

Lu ap (q) =lu + 12tt 2 -^- [1 - exp(-L/12ax)] (13) 

in the size range L > lQxa (fig.l). Of course, when using eq. (fl~3|) in practical 
cases, we should make sure that the quadratic approximation holds down to the 
size dealt with. 

Now, as an example, let us consider the case of Raman shift in titanate 
(TiO-i) nanocrystals for which was proposed the average phonon curve [12, 13] 

u(q) = ojq + A [1 - cos(qa)] (14) 

where ujq = 144 cmT x , A = 20 cm" 1 , and a = 0.376 nm, the linewidth 
being To = 7 cmT 1 so that x = 1.7. Fig. 2 shows the deviations of Raman 
shifts from luq as calculated from eq. ([2|) with ui(q) given in eq. (|14|) (triangles), 
as calculated from eq. ((8|) (circles) and as calculated from eq. (|13|) (dashed 
curve). As expected from Fig.l, in the range of small sizes, the latter equations 
underestimates the deviations of shift. However, departure from the exact curve 
remains relatively small down to 5 nm where the PC model is believed to become 
unreliable [6,12,13]. 

A further example is given by the case of silicon nanocrystals whose phonon 
curve, established by neutron scattering measurements [14], is 

oj{q) = [B 1 +B 2 cos{qa/4)} 1/2 

with a = 0.5483nm, B\ = 1.714 10 5 cmT 2 , B2 = 10 5 cto~ 2 , corresponding 
to ujq = 520.5 cmT 1 and, in the quadratic approximation, A — —6 cmT x \ the 
linewidth is Tq = 3.6 cmT 1 so that x = 1.3. In this case deviations have a 
reversed sign. Departures from the exact curve have roughly the same relative 
values found in the previous case. It can be noted that because of the factor 
1/4 in the cosine argument, the quadratic approximation appears to hold in the 
whole range considered. 

As a final consideration we note from Fig. 1 that one could improve the 
approximation of curve (jT3)) (cases x> 1) by introducing a suitable factor to 
©,5(77). In this connection, we found that —In [&s(r))] = 0.3?7exp(— rj/2) + n/12 
improves the agreement down to n ~ 6. However, in our opinion, this correction 
does not recompense the loss of simplicity in eq. (|13() . 

In conclusion, we have shown that in the framework of the PC effect as 
described by the hard confinement model, the Raman shift in nanocrystals can 
be calculated, with good approximation, by means of a simple analytical form in 
the size range where the quadratic approximation of the average phonon curve 
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holds. The goal is that the dependence of shifts on model parameters as well as 
on particle sizes is immediately recognizable. 
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Captions of Figures 

Fig.l Negative logarithm of the function 0s(?y)(eq. 10). Inset, dependence 
of — ln[®s{v)] on x f° r small r\: a, x=0.5; b, x=l; c, x=2; d, x=3. Dashed 
curves represent parabola with equations / (77) = ■mx 2 T] 2 : a, m=0.6; b, m=0.8; 
c, m=l; d, m=l. 



Fig. 2. Deviations of Raman shifts in nanocrystals from the Raman shift 
in the infinite crystal. Curves are referred to the source equations. The upper 
curves have been calculated for the Titanatc and the lower curves for the Silicon. 
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